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A matrix Hamiltonian model is developed to address the finite-volume effects appearing in studies of baryon 
resonances in lattice QCD. The Hamiltonian model includes interaction terms in a transparent way and can be 
readily generalized to address multichannel problems. The eigenvalue equation of the model is exactly solvable 
and can be matched onto chiral effective field theory. The model is investigated in the case of A — > Nn 
scattering. A robust method for determining the resonance parameters from lattice QCD is developed. It involves 
constraining the free parameters of the model based on the lattice spectrum in question. The method is tested 
in the context of a set of pseudodata, and a picture of the model dependence is obtained by examining a variety 
of regularization schemes in the model. A comparison is made with the Liischer method, and it is found that 
the matrix Hamiltonian method is equally robust. Both methods are tested in a more realistic scenario, where a 
background interaction corresponding to direct Nir 4-> Nn scattering is incorporated into the pseudodata. The 
resulting extraction of the resonance parameters associated with the A baryon resonance provides evidence that 
an effective field theory style of approach yields a successful realization of finite- volume effects in the context 
of baryon resonances. 

PACS numbers: 12.38.Gc 12.38.Aw 12.39.Fe 



I. INTRODUCTION 

The behaviour of the excited states of baryons represents an 
important research topic in nuclear and particle physics. The 
extraction of resonance parameters (e.g. masses and widths) 
provides insight into the scattering behaviour of hadronic in- 
teractions. Lattice QCD constitutes a principal tool for the 
analysis of nonperturbative physics; however, observables 
must be computed in a finite volume. Finite- volume models 
of hadronic interactions represent a valuable avenue for inves- 
tigating the intrinsic features of lattice QCD. They assist in the 
interpretation of discrete eigenstates in relation to the asymp- 
totic behaviour of the S matrix measured in experiments. Re- 
cent advances in lattice QCD have enabled the extraction of 
the lowest-lying excitation energies [1-5]. However, the ex- 
traction of resonance parameters from multihadron interac- 
tions presents an ongoing challenge for current research [6- 
16]. 

In this article, a finite-volume matrix Hamiltonian model is 
introduced, which is exactly solvable. A matrix is constructed, 
and its entries are populated with energies corresponding 
to different momentum states of the interaction channel(s). 
The solution of the eigenvalue equation takes the form of a 
one-loop renormalization formula for a bare resonance en- 
ergy, reminiscent of finite-volume chiral effective field the- 
ory (%EFT). The Hamiltonian can then be matched onto the 
corresponding xEFT formula by choosing the form of the in- 
teraction in the model. A property of the formalism is the abil- 
ity to reproduce the "avoided level crossing" observed in lat- 
tice QCD calculations, while also being able to reproduce the 
appropriate continuum limit of xEFT. Thus the physical be- 
haviour of the energy eigenvalues near a resonance due to the 
mixing of the particle states can be estimated. Furthermore, 
the Hamiltonian method lends itself to an intuitive generaliza- 
tion in addressing multichannel problems, with the inclusion 



of additional degrees of freedom, corresponding to the new 
interactions. 

The matrix Hamiltonian formalism is applied to A — > Nir 
decay, and its energy spectrum is generated for a range of 
finite box sizes, L. A method for extracting the resonance 
parameters from a lattice QCD spectrum is then postulated, 
which involves fitting the free parameters of the model. The 
method is tested by generating a set of pseudodata from the 
model itself. The pseudodata are analysed with an alternative 
version of the model and a selection of different regulariza- 
tion parameters. This provides an indication of the model de- 
pendence in correctly obtaining the resonance position. The 
success of the identification of the phase shift is compared to 
that of the Liischer method [17]. This serves to establish a 
benchmark for the Hamiltonian model in obtaining the reso- 
nance position. It is found that the matrix Hamiltonian model, 
as applied to A — > Nit decay, is of comparable accuracy to 
the Liischer method. 



H. SCATTERING THEORY FROM EFFECTIVE FIELD 
THEORY 



In continuum scattering theory, below inelastic thresholds, 
amplitudes of the partial-wave decomposition of the wave- 
function, which are complex valued, may be converted into 
a real-valued number called the phase shift, 5i (k), by exploit- 
ing the conservation of angular momentum. The phase shift, 
which depends on the orbital angular momentum I and ex- 
ternal momentum k, is related to the total cross section of a 
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scattering process, via 

a(k) = Jdn\f(k,9)\ 2 , 
where f(k, 6) = ^{21 + l)Pi(cos6) 



(1) 



e iSi(k) S i n ^(fc) 
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where Pi is the Zth Legendre polynomial. The external mo- 
mentum of a particle with mass m is related to the external 
energy by: k — y/E 2 — to 2 . In the case of P-wave scattering 
(I = 1), the total cross section simplifies to the form 



cr(fc) = 47T 



3sin 2 ( 5(fc) 
k 2 ; 



(3) 



where S(k) = Si(k). The energy at which the phase shift 
passes through 90° is the resonance energy, E ICS . 

The investigation of A — > Ntt decay is of particular in- 
terest. It provides insight into the A resonance of Ntt scat- 
tering; an archetypal example of a baryon resonance. The on- 
shell t matrix associated with Ntt scattering from a A baryon, 
T = t(k,k; E + ), can be obtained by solving the Lippmann- 
Schwinger equation, and the phase shift is directly related to 
the t matrix, via 
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(4) 
(5) 



The quantity A represents the bare value of the resonance 
energy, which becomes renormalized by the one-pion loop in- 
tegral, Satv(^) (shown in Fig. 1), and may be chosen such 
that _E rcs matches the phenomenological value, Ef^/ S = 292 
MeV. By working in the heavy-baryon approximation, the ki- 
netic energy of the nucleon is neglected in this simple model, 
and calculations may be performed relative to the nucleon 
mass. Thus, the external energy may be expressed in terms 
of the pion energy, E = w w (fc). 

The coupling gAN{k) takes the following form, as obtained 
from xEFT, 



ffAivO) = XA 



2 k 2 u 2 {k) 
TT u n (h) 



(6) 



where the coefficient xa is defined in terms of known phe- 
nomenological parameters: f n = 92.4 MeV and C = —1.52, 
as derived from the SU(6) flavour-symmetry relation 



XA = 
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(7) 



The parameters / w and C occur in the first-order interaction 
Lagrangian of chiral perturbation theory (xPT) describing the 
A — > Ntt interaction [18-23] 



: x pt = -—NT a 7 ^A v d„w a + h.c, (8) 
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FIG. 1: The leading-order one-pion loop contribution to the reso- 
nance energy E associated with A — > Ntt decay. All charge con- 
serving transitions are implicit. 



where T a are the relevant 2x4 isospin transition matrices 
[24]. 

Note that a finite-range regulator u(k) has been introduced 
into the coupling term of Eq. (6). This serves to control the 
otherwise ultraviolet divergent loop integral and, in addition, 
will be used in Sec. Ill to establish a finite interaction range 
in the Hamiltonian model. For more detailed discussions of 
finite-range regularization, see Refs. [25-29]. 

The leading-order pion loop integral shown in Fig. 1 takes 
the form 
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where V indicates that a principal value integral must be 
performed. By using Eq. (5), it is straightforward to solve 
for (5(fc). A plot of the phase shift against external energy, 
E = is shown in Fig. 2, which takes the form of a 

Breit-Wigner-like curve. 



III. FINITE- VOLUME HAMILTONIAN MODEL 

In the following finite-volume heavy-baryon model, a ma- 
trix Hamiltonian is defined to include the interaction between 
the A baryon and the pion-nucleon system. The rows and 
columns of H represent the total three-momentum values 
available to the pion (fci = 2w/L, fc 2 = 47r/L . . .). In the 
simplest version of the model, the Ntt states only couple to 
the A baryon. 

The Hamiltonian may be written as the sum of free and in- 
teraction Hamiltonians, H = Ho+Hj. The free Hamiltonian, 
Ho, takes the form 



Ho 



/A 
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where w,({;„) = \/k 2 + m 2 is the pion energy, and the 
pion momenta at finite volume take only the discrete val- 
ues available in a box of length L, obeying the condition: 
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FIG. 2: (color online). The infinite-volume phase shift 5 associated 
with elastic Nn scattering via a A baryon intermediate state, plotted 
against the external pion energy E (where E to t ~ Mjv + E), as 
calculated from the on-shell t matrix of Eq. (5). 



kl = ) 



n, for a squared integer n 



has nonzero diagonal entries representing the noninteracting 
energy of a pion at certain values of momentum. The nucleon 
recoil energy vanishes in the heavy-baryon limit. The first ele- 
ment of Hq represents the bare resonance mass relative to the 
nucleon mass, Ao. The value of Ao is chosen so that the res- 
onance position takes the phenomenological value 292 MeV, 
in the infinite-volume case in Eq. (4). 

The interacting Hamiltonian Hj contains couplings g&N in 
the top row and left-most column, 



Hj 
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(12) 



This represents the interaction between the bare A baryon and 
the Ntt system over a range of momentum values. The cou- 
plings are chosen to take a similar form to that of Eq. (6) but 
contain additional factors relevant to the normalization of the 
discrete momenta in a finite-volume box, 



9 AN (k n ) 
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Note that the cubic symmetry of the system means that each 
nonzero n need occur only once in the rows and columns of 
H, so long as the finite-volume couplings g% u N {k n ) contain 
the appropriate normalization factor \/ Cj,(n), where C^ip) 
represents the number of ways of summing the squares of 
three integers to equal n. 

For the purposes of this model, u{k n ) is chosen to take the 
form of a dipole regulator, with a regularization scale of A = 
0.8 GeV, 



u(k n ) 



k 2 
A 2 



(14) 



However, in Sec. V, the associated model dependence of the 
choice of regulator on the extraction of the resonance param- 
eters from the model is investigated and a variety of regular- 
ization schemes are examined. 

The eigenvalues of the Hamiltonian matrix correspond to 
the different energy levels in a finite volume, and may be cal- 
culated for a range of L values. The ten lowest energy levels 
as a function of L are shown in Fig. 3, with the noninteracting 
energies shown as dotted lines. 

The Hamiltonian model, by construction, matches onto the 
XEFT formula for the renormalization of the resonance mass 
in a finite volume. The eigenvalue equation of the Hamilto- 
nian, det(H — XT) = 0, can be solved exactly, and takes the 
form 



A = Ar 



Z -' cJ„-(fc„) - A 
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A comparison with Eq. (10) shows that the sum term takes 
the same form as the loop integral, upon the restoration of the 
continuum integral 



-r >Mn) 



n=l 



d s k. 



(17) 



This correspondence is a consequence of the judicious choice 
of the form of the coupling in the finite-volume Hamiltonian 

_fin 
"AiV' 

The subtle, yet important difference between Eqs. (10) and 
(16) is that the energy eigenvalue A appears on both sides of 
the equation derived from the Hamiltonian. The presence of 
A in the denominator, whose solution must be finite, results 
in the avoided level crossing observed in lattice QCD, where 
the resonance energy corresponds to maximal mixing of the 
single-particle and multiparticle states. That is, the correct 
quantum mechanical behaviour of the multiparticle system is 
directly built into the Hamiltonian model. This is the finite- 
volume analogue of the principal value integral in the contin- 
uum theory. 

The Hamiltonian model may be extended to include a direct 
Ntt O Ntt interaction, which could be viewed as an approx- 
imation to the Chew-Low interaction. For example, one can 
add additional terms in the interaction matrix, Hj: 



gcL(ki,kj) = -gN-n{h)gNTi{kj) 



1 



(18) 



In the cloudy bag model [30, 31], the couplings g^^ and g^ N 
are related by the equation gNTr(ki) — Kg^ N (ki) with a ratio 
k = -y/25/18 estimated from SU(6) symmetry. In the case of 
the t matrix, the Lippmann-Schwinger equation with a rank- 2 
separable potential can be solved as per Mongan [32]. 

This additional interaction will be used to test the model 
independence in the procedure of extracting the phase shift in 
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FIG. 3: (color online). The lowest-lying energy levels from the ANn 
model, described in Eqs. (1 1)— (14) (solid lines), and the correspond- 
ing noninteracting energies (dotted lines). 



Sec. V. The comparison of two models, one with extra back- 
ground interactions of the form of Eq. (18) and one without 
the extra interactions, will provide insight into the robustness 
of the identification of resonance parameters from an effective 
model. 



IV. PHASE SHIFT EXTRACTION FROM THE LATTICE 

The identification of resonance parameters in a finite vol- 
ume is important for interpreting lattice QCD results in the 
light of experiment. The connection between the energy spec- 
tra of lattice QCD simulations and the asymptotic states of the 
S matrix represents a current challenge for research in finite - 
volume hadron spectroscopy [6, 33-36]. 

The features of the Hamiltonian model are ideally suited to 
the extraction of resonance parameters from lattice QCD en- 
ergy spectra. Not only is the construction of the Hamiltonian 
transparent and intuitive in terms of the potentials one chooses 
to include, but the formalism may be easily generalized to in- 
clude multiple channels. By generating energy eigenvalues 
from a Hamiltonian model, such as those shown in Fig. 3 for 
the A — ► Ntt interaction, the free parameters of the model 
may be tuned to fit the energy spectrum of a lattice QCD cal- 
culation. In the case of the A — > Ntt model, this entails min- 
imizing the chi-square for the parameters \A and Ao. The 
resonance energy and phase shift may then be extracted by us- 
ing the fitted values of the parameters, input into Eq. (5). An 
infinite-volume phase shift, much like that of Fig. 2, is thus 
obtained, except that the underlying parameters have been ex- 
tracted directly from lattice QCD results. 



V. COMPARISON WITH THE LUSCHER METHOD 

In order to demonstrate the success of the Hamiltonian 
method in reproducing accurate phase shifts from the finite- 
volume energy spectrum of the model, a comparison is made 



with the well-known Liischer method [17]. Liischer's formula 
describes a fixed relationship between the scattering phase 
shift S and the momentum corresponding to the jth energy 
level in a finite volume 



S(kj-,L) = jir-(f> 



kjL 
2tt 



(19) 



The kinematic function <fi(q) takes the form of a three- 
dimensional Zeta-like function (which must be regularized in 
an appropriate fashion), defined in terms of dimensionless lat- 
tice momenta q= kL/ (2ir), 



4>(q) = — arctan 
Zoo(l,q 2 ) 



(20) 
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The momenta corresponding to a lattice energy spectrum may 
be input into Liischer's formula to obtain volume-dependent 
"phase shifts" 5(kj; L). 

As an example, consider the application of the Liischer 
method to the energy spectrum of the A — > Ntt model, treat- 
ing the spectrum as pseudodata. By taking the momenta asso- 
ciated with the energy eigenvalues, shown in Fig. 3, as inputs, 
and choosing a particular box size L, the finite-volume esti- 
mates of the phase shifts may be compared with the known 
infinite-volume phase shifts, as shown in Fig. 2. The result 
is plotted on the same axes, in Fig. 4, for a range of box 
sizes. Both the infinite-volume phase shift and the Hamil- 
tonian energy spectrum are generated using a dipole regula- 
tor with A = 0.8 GeV. Figure 4 shows that the finite-volume 
corrections associated with applying the Liischer method to 
a finite-volume effective Hamiltonian model are insignificant 
for L > 3 fm at the physical pion mass. It will be demon- 
strated that the matrix Hamiltonian method of obtaining a 
phase shift from finite-volume energy levels is comparable 
with the Liischer method. 

Note that, in using Liischer's method, an interpolation func- 
tion must be chosen in order to obtain the pole position from 
the phase shift. The Breit-Wigner curve represents a suitable 
functional form, 



l5{k) sm5(k) 
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E - E K 
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(22) 



where the width contains a cubic momentum factor, which 
yields the correct phase-space dependence for I = 1 [37], 



£.3 



(23) 



To illustrate the comparison between Liischer's method and 
the matrix Hamiltonian method, the pseudodata obtained from 
the A — !• Ntt model will again be used. A modified version of 
the A — > Ntt model is constructed by replacing the regulator 
function u(k n ) with a different form of regularization. Firstly, 
a Gaussian regulator with a variety of values for the regular- 
ization scale A will be considered. Using the A — > Ntt model, 
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FIG. 4: (color online). Estimates of the phase shift obtained from Liischer's 
formula applied to the A — > Nn model for a range of L values. The infinite- 
volume phase shift is also shown as a solid line. 



the closest two energy eigenvalues to the value of E res , as es- 
timated by Liischer's formula, are chosen in order to constrain 
XA and An. 

The behaviour of the resonance energy E lcs may be plotted 
as a function of 1/L, as shown in Fig. 5. A key observation 
is that the periodic deviations in the extraction of E rcs arise 
from the fact that, at certain volumes, an energy level may or 
may not lie near the resonance position (where 8 = 90°). That 
is, they are simply artefacts associated with the choice of in- 
terpolation. For box sizes where a discrete energy eigenvalue 
does lie near E rcs , the experimental value is attained within a 
few MeV. This can provide a guide in choosing suitable box 
sizes for reliably reproducing the position of the A resonance. 

It is also worth noting that the use of a Gaussian regula- 
tor, with A ~ 0.6 GeV, provides a much more reliable re- 
construction of the resonance position from the pseudodata 
than a low or high value of A. This is a consequence of a 
close match between the underlying regulator of the pseudo- 
data, and the Gaussian regulator at this value of A. It provides 
evidence that a choice of regulator that correctly reflects the 
underlying phenomenology leads to a more stable identifica- 
tion of the resonance position. Nevertheless, even in the case 
of A — > oo, which essentially corresponds to a dimensional 
regularization-like result [26], the matrix Hamiltonian method 
is just as successful as the Luscher method. 

Note that in Fig. 4 use of the Luscher method results in 
an excellent agreement with the infinite-volume phase shift. 
However, at small volumes, the sparse distribution of energy 
eigenvalues leads to a difficulty in extracting the behaviour at 
8 = 90°. This accounts for the model dependence apparent in 
Fig. 5. 

One might argue that the analysis of the pseudodata using 
a similar model that includes only the A — > Nir couplings is 
ideally suited to the reconstruction method proposed above. 
In a more realistic scenario, a new pseudodata set is generated 
from a Hamiltonian model that now includes the effective con- 
tact interaction, introduced in Eq. (18). Again, a dipole regu- 
lator with A = 0.8 GeV is used. The energy levels of the new 
pseudodata are fitted to the more basic model using a Gaus- 
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FIG. 5: (color online). The resonance energy E rcs plotted against 
1/L. The experimental value is marked with a square. The results 
from fitting pseudodata with a different regulator, a Gaussian with 
A = 0.5, 0.6 and 0.8 GeV, or treating A as a fit parameter, are plot- 
ted. The result of effectively removing the regulator (i.e. A — 5* oo) is 
also shown. For comparison, the approach using Liischer's method 
is marked with a solid line. 
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FIG. 6: (color online). The resonance energy, E rBS , plotted against 
1/L, resulting from fitting pseudodata containing the nonresonant 
background interaction of Eq. (18). The experimental value is 
marked with a square. The result from fitting pseudodata with a 
Gaussian regulator is plotted, treating A as a fit parameter. The re- 
sult corresponding to A — > oo is also shown. Liischer's approach is 
marked with a solid line. 



sian regulator, and the regularization scale, A, is treated as an 
additional fit parameter. The behaviour of the resonance en- 
ergy versus 1/L is displayed in Fig. 6. The fitted values of A 
in Fig. 5 are recalculated at each value of L, but they are typi- 
cally close to 0.6 GeV, whereas in Fig. 6 the values are closer 
to 0.5 GeV. The tendency for slightly smaller fitted values of 
A in Fig. 6 simply reflects the compensation of the regulator 
for the missing background interaction. For all regularization 
methods tested, the presence of a background interaction also 
causes an exaggeration in the periodic deviations of E rcs from 
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interpolation. In all cases, the deviations are much larger due 
to the effect of the nonresonant background interaction. Nev- 
ertheless, a fairly reliable extraction of E res is possible, if the 
energy eigenvalues lie close enough to the resonance energy. 
Most notably, there is significant improvement in the identi- 
fication of E les using the matrix Hamiltonian method if A is 
treated as a fit parameter. This is remarkable, as it shows that 
the bare resonance mass, the coefficient of the interaction and 
the finite-range scale A are collectively able to compensate for 
the missing Ntt O Ntt interactions; the underlying interac- 
tions are simply incorporated into the fit degrees of freedom. 



VI. ANALYSIS OF THE EIGENVECTORS 

An analysis of the eigenvectors serves to emphasise a key 
feature of the Hamiltonian model. While only a single local 
operator describing the A resonance is included in the model, 
the interactions drive a significant coupling to all nearby lat- 
tice eigenstates. That is, one observes the generation of quan- 
tum mechanical admixtures of states from the interactions. 

Consider the square of the overlap of each eigenvector with 
the bare state | Ao), with a volume-dependent coefficient cal- 
culated from the following Jacobian: 
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The result may be plotted as a function of external energy E. 
The peak of the resulting curve should lie at the value of the 
renormalized mass of the A baryon relative to the nucleon, 
292 MeV. The result for the matrix Hamiltonian model, using 
a dipole regulator with A = 0.8 Ge V, and including the Chew- 
Low-like background interactions from Eq. (18), is shown in 
Fig. 7. The magnitude of the contribution from the eigenvec- 
tor spectrum to Ao, |(Ao|v)| 2 , is shown for box sizes of 4, 8 
and 16 fm. Over a broad range of energies, which dominate 
across the width of the resonance, it is evident that no single 
eigenstate, by itself, can be identified as either a local or a mul- 
tihadron state. That is, the local, noninteracting Ao state cou- 
ples to a range of low-lying states. This illustrates how early 
ideas from the heavier quark-mass domain, in distinguishing 
single-particle and multiparticle states, are no longer applica- 
ble in the light-quark regime of modern lattice simulations. 

VII. CONCLUSION 

Finite-volume models provide key insights into the intrinsic 
attributes of lattice QCD. By examining the behaviour of res- 
onances and multiparticle states, one is able to identify more 
precisely the impact on a calculation associated with finite- 
volume effects. 

By introducing a matrix Hamiltonian method, the finite- 
volume behaviour near a resonance is reproducible. The 
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FIG. 7: (color online). The contribution of the eigenvectors to the bare A 
state in a finite volume, for a range of box sizes. The finite-volume Hamilto- 
nian used includes the Nn +4- Nn interaction terms. 



avoided level crossings that indicate nontrivial mixing be- 
tween single-particle and multiparticle states are also present 
in the model. Because the Hamiltonian is exactly solvable, 
the eigenvalue equation may be matched onto chiral effective 
field theory simply by a choice of coupling parameter for the 
potential. The Hamiltonian model can also be easily gener- 
alized to include multiple channels, in order to address more 
difficult scattering problems. 

As a simple, introductory example, the A baryon resonance 
of Ntt scattering was considered, and a method for extracting 
the phase shift from lattice QCD was developed. By matching 
the energy eigenstates of the model to those of a lattice cal- 
culation, the free parameters of the model were constrained, 
and their values were input into a scattering t matrix equa- 
tion, which uses chiral effective field theory. Initially, only 
the leading-order one-pion loop integral contribution to the 
resonance energy was included. 

The method was tested in the context of a set of pseudodata 
generated by the A — > Ntt model. By considering a variety 
of different regularization scales, an indication of the model 
dependence in the calculation was found. A comparison was 
made with the well-known Liischer method, and it was discov- 
ered that the matrix Hamiltonian method was at least compa- 
rable in reliability, with improvement when the regulator is a 
good approximation of the underlying phenomenology. 

The extraction of the resonance parameters was investi- 
gated in a more realistic context, where nonresonant back- 
ground interaction terms (corresponding to direct Ntt ■<-» Ntt 
scattering) were incorporated into the pseudodata. It was 
found that the A — > Ntt model was able to improve upon 
the Liischer method, if the regularization scale was treated as 
an additional fit parameter. This indicates the ability of a bare 
resonance mass, a coupling scale, and regulator to compen- 
sate for underlying interactions that do not appear directly in 
the model. It also indicates that more complicated interactions 
may be effectively incorporated into the degrees of freedom of 
the model. 

Thus the matrix Hamiltonian model, which uses an ef- 
fective field theory style of approach, may indeed act as a 
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guide to lattice QCD calculations in identifying the salient 
features derived from finite-volume effects, and the impact 
of the leading-order chiral contributions near a resonance. 
Though lattice QCD calculations are typically performed at 
pion masses heavier than the physical value, the avoided level 
crossing feature is expected to occur on the lattice when the 
pion-nucleon threshold appears below the A resonance. After 
considering some ideal pseudodata, the matrix Hamiltonian 
approach is well equipped to describe the behaviour of actual 
lattice QCD results. 
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